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1.  Summary  of  Research  on  Grant  AFOSR-79-0029 

Research  supported  by  this  grant  has  been  concentrated  in  four  distinct,  but 
related  projects:  1)  an  analysis  of  inherent  phase  distortion  in  rectangular 
and  polar  format  FFT  processing  algorithms,  2)  a  SAR  computer  simulation  with 
polar  format  recording,  3)  a  study  of  2D  interpolators  for  polar-to-rectangular 
coordinate  transformation,  and  4)  an  investigation  of  number  theoretic  concepts 
for  high  speed  failure  resistant  digital  processors  required  in  real-time  SAR 
systems.  The  following  paragraphs  summarize  results  in  these  four  areas. 

1.1  Analysis  of  Inherent  Phase  Distortion 

Our  study  of  SAR  signal  processing  algorithms  began  with  an  analytical 
evaluation  of  inherent  phase  errors  that  occur  in  "stretclv'  processing  that  is 
popular  in  present  day  spotlight  mode  SARS.  In  stretch  systems,  linear  FM 
waveforms  are  transmitted  at  N  distinct  positions  along  the  flight  path,  at 
which  points  the  returns  are  received,  demodulated,  sampled,  and  recorded  for 
subsequent  processing.  This  creates  a  linear  FM  variation  in  each  return  from  a 
point  target  along  the  range  coordinate  (time).  The  relative  motion  of  the 
radar  and  the  ground  creates  an  approximate  linear  FM  variation  in  point  target 
responses  along  the  azimuthal  coordinate  (direction  of  flight) ,  due  to  the 
Doppler  effect.  After  the  data  have  been  recorded  at  all  N  positions,  a 
reference  function  is  mixed  with  the  data  to  focus  a  predetermined  reference 
target,  and  the  ensemble  is  processed  witli  a  2D  FFT  to  resolve  the  targets  by 
spectral  analysis.  Stretch  processing  is  motivated  by  the  computational 
efficiency  of  the  FFT,  which  holds  potential  for  real-time  processing. 

Unfortunately,  there  are  two  serious  deficiences  in  stretch  signal  pro¬ 
cessing.  The  first  is  that  the  azimuthal  frequency  variation  is  only  approx¬ 
imately  linear  FM,  so  that  there  are  residual  phase  errors  which  depend  on 
target  position,  and  which  cannot  be  corrected  by  the  FFT  algorithm.  The 
second  is  that  targets  move  in  range  relative  to  the  radar,  so  that  polar  format 
recording  is  required  for  spotlight  mode  operation.  Since  the  FFT  requires 
data  samples  in  a  rectangular  coordinate  system,  2D  interpolation  is  needed  to 
generate  data  samples  which  lie  on  a  rectangular  raster.  The  2D  interpolation 
requires  considerable  computation,  to  the  extent  that  it  may  be  a  dominant 
factor  in  the  overall  computational  complexity  of  the  stretch  algorithm. 


In  our  study  of  phase  error  effects,  analytical  expressions  were  derived 
for  the  phase  of  the  recorded  response  from  point  targets  in  the  ground  plane. 
Phase  functions  were  derived  for  both  polar  and  rectangular  recording  systems. 
Since  the  phase  response  is  a  relatively  complicated  nonlinear  function  of 
target  position,  it  was  expanded  in  a  Taylor  series  so  the  effects  of  linear 
phase  error,  quadratic  phase  error,  and  higher  order  phase  error  terms  could  be 
identified.  The  results  of  this  analysis  isolate  the  effects  of  linear  phase 
error  (misregistration)  and  quadratic  phase  error  (defocus)  in  both  the  polar 
and  rectangular  domains  [20]. 

In  an  independent  analysis,  residual  phase  errors  were  analyzed  for  an 
algorithm  which  uses  FFT  stretch  processing  in  the  azimuthal  coordinate  only, 
assuming  that  range  processing  is  implemented  independently,  possibly  by  real¬ 
time  analog  devices.  The  results  of  this  analysis  are  useful  for  relating  phase 
error  to  image  size  for  systems  using  hybrid  signal  processing  concepts  [17]. 

1.2  SAR  Computer  Simulations 

A  SAR  computer  simulation  was  developed  that  processes  recorded  data  files 
to  simulate  system  response  to  ideal  point  targets.  Artificial  data  files  are 
generated  from  the  analytical  phase  functions  derived  in  project  1.  The  data 
files  are  then  used  as  input  to  a  simulated  SAR  processor  for  either  rectangular 
format  or  polar  format  signal  processing.  The  simulation  is  useful  for  studying 
the  effects  of  residual  phase  errors  and  for  comparing  the  performance  of  differ¬ 
ent  processing  algorithms. 

For  example.  Figure  1  shows  the  response  of  a  spotlight  mode  SAR  for  two 
point  targets  processed  with  rectangular  format  recording.  The  center  target 
is  located  at  the  reference  position,  and  the  response  shown  is  the  ideal  response 
of  a  Kaiser  window.  The  image  of  the  second  target  is  smeared  (defocussed)  and 
misplaced  (misregistered) .  This  example  illustrates  the  serious  degradation  that 
results  from  using  rectangular  format  porcessing  in  spotlight  mode  operation. 
Figure  2  shows  the  response  that  is  obtained  when  ideal  polar  format  recording 
is  used.  The  reference  target  is  again  perfectly  focussed,  while  the  second 
target  is  only  slightly  degraded  by  residual  phase  errors.  Figure  3  shows  the 
system  response  when  the  polar-to-rectangular  coordinate  transformation  is  im¬ 
plemented  with  a  2!)  inverse  distance  (nonseparable)  interpolation  algorithm. 

This  simulation  shows  that  the  smeared  target  can  he  recovered  through  correct 
processing,  although  there  is  evidence  in  Figure  3  of  serious  noise  created 


oint  target  response  with  rectangular  format  record 


point  target  response  with  ideal  polar  format  recording 
tion  error). 
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Figure  3  Simulation  of  point  target  response  with  polar  format  recording  and  1st  order 
nonscparnhle  int orpolat  i  on . 
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by  the  crude  nature  of  the  inverse  interpolation  algorithm.  Simulation  studies 
have  investigated  relationships  between  spread  (defocus)  and  both  target  distance 
from  the  reference,  and  resolution  [20], 

1.3  Two-Dimenst ional  Interpolators 

In  digital  SAR  processing  2D  interpolation  for  polar-to-rectangular  coordin¬ 
ate  conversion  is  very  complicated  because  the  initial  data  points  are  non- 
uniformlv  spaced  relative  to  the  rectangular  coordinate  axes.  Five  types  of  2D 
interpolators  are  presently  under  study:  1)  1st  order  inverse  distance,  2)  1st 
order  separable  bilinear,  3)  3rd  order  separable  Lagrange,  4)  separable  optimal 
FIR  filtering,  and  5)  zero  order  interpolation  with  over-sampling.  These  simple 
interpolators  have  been  selected  because  they  have  low  computational  complexities, 
although  there  are  important  questions  to  be  answered  regarding  their  quality  of 
performance.  Attempts  are  now  in  progress  to  characterize  performance  in  terms  of 
interpolation  error.  For  example.  Figures  4  and  5  show  the  frequency  responses 
of  the  1st  order  inverse  distance  and  separable  bilinear  interpolators,  respec¬ 
tively,  as  applied  to  uniformly  spaced  data  samples.  The  energy  in  the  sidelobes 
is  lower  for  the  separable  bilinear  filter,  implying  that  this  filter  should  re¬ 
sult  in  better  spacial  domain  performance.  An  important  thrust  of  current  work 
is  to  characterize  performance  through  analytical  techniques  which  represent  the 
quality  of  the  interpolator  in  terms  of  the  sidelobe  energy  in  the  spectrum  [19]. 

1.4  Residue  Number  Concepts  for  SAR  Digital  Processors 

Research  on  this  project  involves  an  investigation  of  residue  number  concepts 
for  realizing  SAR  digital  processors  with  high  speed  capability  and  failure  re¬ 
sistant  behavior.  This  work  is  motivated  because  real-time  SAR  processing  is 
severely  demanding  on  the  resources  of  a  digital  system,  to  the  extent  that 
dedicated  processors  are  needed  to  satisfy  data  rate  requirements.  Residue 
number  artithmetic  lias  been  known  for  many  years  for  its  potential  for  high 
speed  addition  and  multiplication,  although  it  has  not  been  effective  in  gen¬ 
eral  purpose  computers  because  some  functions,  such  as  sign  detection,  magni¬ 
tude  comparison,  and  division  are  awkward  and  time  consuming  operators  in  the 
residue  code.  Within  the  last  few  years  residue  number  arithmetic  has  captured 
the  attention  of  researchers  who  are  interested  in  digital  signal  processing 
(digital  filtering,  spectal  analysis,  docimution/interpolation,  statistical 
analysis,  etc.)  rather  than  in  the  broader  area  of  general  purpose  computing 
111,  151. 
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In  addition  to  high  speed  capabilities,  residue  number  codes  have  a  mathe¬ 
matical  structure  that  facilitates  error  detection  and  correction  when  proper 
recundancy  is  added.  These  properties  result  from  the  parallel  residue  structure, 
(which  prevents  the  propagation  of  errors)  and  the  lack  of  a  hierarchical  signif¬ 
icance  among  the  digits, (which  allows  erroneous  digits  to  be  discarded  without 
rendering  the  result  completely  useless).  Our  recent  work  at  the  Coordinated 
Science  Laboratory  has  examined  the  theory  of  systematic  redundant  RNS  coding 
for  the  isolation  and  detection  of  hardware  failures,  with  subsequent  "soft  failure" 
that  allows  a  faulty  processor  to  continue  operating  with  reduced  capabilities 
[5,  7,  8,  14].  The  theory  is  elegant  and  its  capabilities  are  attractive  in 
relation  to  the  special  requirements  of  SAR  signal  processing.  One  major  concern 
at  this  time  involves  hardware  complexity  in  the  error  checkers  and  the  persis¬ 
tent  question  of  reliability  in  the  error  checking  hardware.  One  of  the  important 
thrusts  of  our  present  work  is  to  determine  desirable  RNS  sturctures  and  im¬ 
plementation  techniques  that  can  be  realized  in  VLSI  technology  [5], 

2.  Future  Directions 

During  our  recent  work  in  SAR  several  general  observations  have  been  made 
which  provide  motivation  and  direction  for  futher  research.  The  first  is  that 
it  has  become  apparent  to  us  that  SAR  processing  is  a  specific  case  of  a  more 
general  class  of  2D  digital  array  processing  algorithms.  The  general  class 
includes  techniques  used  in  computer-aided  tomography,  beam  forming  sonar, 
geophysical  s^emic  processing,  and  radio  astronomy.  It  is  our  belief  that 
future  SAR  studies  should  include  an  evaluation  and  analysis  of  alternative 
2D  array  algorithms  used  in  these  related  areas.  A  second  general  observation 
is  that  many  2D  array  systems  have  one  dimension  that  is  inherently  analog, 
and  a  second  dimension  that  is  inherently  discrete.  It  is  our  belief  that 
future  studies  should  include  hybrid  signal  processing  techniques  that  will 
use  both  modern  analog  devices  and  sophisticated  digital  techniques  to  their 
maximum  capabilities.  Our  third  general  observation  is  that  failure  resistance 
by  means  of  residue  number  coding  may  be  more  efficiently  realized  by  using  the 
class  of  RMS  product  codes.  Virtually  all  present  work  in  redundant  RNS  signal 
processing  has  concentrated  on  the  class  of  systematic  RNS  codes.  There  is  some 
evidence  that  the  product  codes  can  lead  to  considerably  less  complexity  in  the 
error  checkers,  and  therefore  will  be  more  useful  in  VLSI  implementation  of 
failure  resistant  processors. 
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